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ABSTRACT 



The general differential equations describing the behavior of the 
rotary regenerator including longitudinal heat conduction in the direction 
of fluid flow are sufficiently complicated to preclude a complete analyti- 
cal solution. A. fev» solutions found in the literature restricted to very 
special cases are discussed. 

A numerical finite difference method is presented which will 
determine the effect of longitudinal heat conduction in rotary regenerators 
for steady state conditions. In the development no assumptions are made 
which would restrict the range of parameters for which the analysis 
would be applicable. 

The conduction effect on the regenerator effectiveness was evaluated 
with the view of obtaining results most useful for the gas turbine 
regenerator problem, however, these results may also be used for other 
regenerator problems. 

A CDC 1604 digital computer was used to carry out the computations. 

The results are presented graphically and in tabular form, employing a 
suitable set of non dimensional parameters. The range of parameters 
which have been covered arc: 

0-5 Ciain/pBiax =^1 -0 
1.0 , Cr/Cn^lO.O 
1.0 NTUo 4 20.0 
0.25 4 (hA)*j 1.: 

0 \ . 0.2 

0.15 As* 1.3 

The author expresses his appreciation to C. P. Howard, Associate 
Professor, for his direction and encouragement in this work. 
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I. Introduction 

The regenerative type of heat exchanger has been subjected to 
extensive treatment within the last thirty years or so and contributions 
to the theory of regenerators have been made by many writers. The most 
questionable of the idealizations made in the development of the theory 
i3 that of zero matrix conductivity in the direction of fluid flow. The 
loos due to heat conduction in the matrix material is specially important 
in the case of high performance regenerators in which a small reduction 
in the effectiveness will result in a large increase of the heat transfer 
areas to compensate for the loss. Another important factor is the modern 
trend toward smaller sizes, since the loss due to conduction of heat in 
the matrix obviously increases with decreasing length. 

r 1 

Schultz 2 solved the problem including conduction for the 
special case of the balanced regenerator and very large values of Cr/C^n 
(i.e., the capacity rate ratio of rotor matrix to minimum fluid capacity 
rate). In this case the system of partial differential equations defining 
the problem for the steady state are reduced to ordinary differential 
equations and the problem can be solved analytically. 

Hahnemann 3 extended the solution to the unbalanced regenerator 
(i.e., unequal partitions of the cross sections and unequal heat transfer 
coefficients together with unequal capacity rate ratios of the fluid 
streams) but still for an infinitely large value of Cr^^in* 

London has shown , as reported by Lambertson i , that a correction 
factor for conductivity in the direction of fluid flow for approximately 

^■Number in brackets refers to bibliography on page 2.1 
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equal fluid capacities, of the form: 

/L\ tg _ i-t A 

C L 

should result in a somewhat pessimistic prediction of the reduction in 
effectiveness. 

It is the objective of this thesis to develop and present a finite 
difference numerical analysis with solutions to solve for the conduction 
effect on effectiveness for the rotary regenerator. Since there arc six 
dimensionless parameters required to specify behavior in the most general 
case, the ranges of the parameters covered in the solutions will be 
those of most interest to the gas turbine regenerator, however, the 
analysis will place no restrictions on these parameters. 
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2. Method 



The conventional idealizations and boundary conditions assumed in 
the derivation of the governing differential equations, including the 
heat conduction in the direction of fluid flow 4 , are the following: 

(See Fig. la) 

a. The thermal conductivity of the matrix material is zero in the 
direction of matrix metal flow. It has a finite value in the 
direction of fluid flo\/, and is infinite in the other direction 
normal to the fluid flow. 

b. The specific heats of the two fluids and matrix material are 
constant with temperature. 

c. Ho leakage of the fluids occurs either due to direct leakage or 
carry over, and each fluid flow is unmixed. 

d. The convective conductance between the fluid and the matrix are 
constant with flow length. 

e. The fluids pass in counter flow directions. 

f. Entering fluid temperatures are uniform over the flow inlet cross 
section and constant with time. 

g. The matrix temperature gradient in the direction of fluid flow is 
zero at the ends. 

h. Regular periodic conditions are established for all matrix 
elements, i.e., steady state condition. 

The differential equations derived by Schultz \2 are given in 
Appendix 1. 

It was suggested b\ IX’.sVa^rra md carried out by Lambertson 1 for 
the case of no conduction, that a rotary regenerator could be divided 
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Cmax 



Typical 



element 



FIG. la ILLUSTRATIVE MATRIX ARRANGEMENT 



AND FLUID FLOW 



into finite elements as represented schematically in Fig. lb. This same 
method can be extended to the case with conduction, if the elements are 
considered to be fixed ir. space, each element can be regarded as a cross 
flow heat exchanger with a gas stream and a metal stream. For an element 
on the side of the heat transfer rate by convection is equal to the 

rate of change in enthalpy of the fluid across the element, 
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( 1 ) 



where _iT av g represents the mean temperature difference between the fluid 
and the matrix element. For a snail enough element the arithmetic mean 
temperature difference may be assumed valid, so that 
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where the fluid and matrix temperatures are subscripted to indicate the 
average temperatures across the inlet and outlet of the element for 
simple unmixed cross flow. 

Considering an energy balance for the element, the energy transferred 
to the element by convection plus energy transferred to the element by 
conduction less the energy transferred out of the element by conduction 
must equal the energy stored in the element: 
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FIG. lb SCHEMATIC REPRESENTATION OF A ROTARY REGENERATOR 



By using equations (l) and (2) the outlet temperatures for an 

element on the side of may be solved for in terms of the remaining 

temperatures and the result expressed in dimensionless form (Appendix 1) 
Since the longitudinal heat conduction is zero at the ends, the 
energy balance equation is different from equation (2) for the elements 
of the first and last row. Three different expressions for the matrix 
outlet temperatures are obtained: 

First row elements 
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Middle row elements 
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Lsst row elements 




(3) 



(A) 



By solving for the fluid outlet temperature on the side of C max in terms 

of the matrix outlet temperature from equation (1), only one expression 

is obtained, which holds for all the elements on the side of C , 

7 max 
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Energy balance; 
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(3) 

By the same method the outlet temperatures of the elements on the side 
of may be obtained using equations (7) and (8). 

First row elements 



, sKi i 



Middle row elements 
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(10) 
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Last vow elements 



(U) 



Solving for the fluid outlet temperature for the elements on the side 
of from equation (7) gives 



which holds for ail the elements on the side of Cpii n . For the values 
of the constants see Appendix 1. 

In addition to the. four dimensionless groupings for no conduction, 
Cmin/Cmax = capacity rate ratio of fluid flow streams 
Cr/C^in ** capacity rate ratio of rotor matrix to minimum 
fluid capacity rate 



two additional dimensionless groupings are needed in the coefficients 
of equations (3), (4), (5), (6), (9), (10), (11), and (12) to take 
conduction into account. These conduction parameters may be defined as: 



since the thermal conductivity and the length of the matrix are the same 
for both sides, and where 




( 12 ) 



•J|f 

(hA) = conductance ratio 



NTUo = 



over-all number of transfer units; 
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conduction parameter on side of C 



'max 






conduction parameter on side of 



Also 



\ 



c 



_ — — — , conduction area ratio. 

A x 

Another form of defining the conduction parameter is 



_ , total conduction parameter 

V. ( , 

In this case As is the total solid area available for conduction, that is 




This form offers the advantage that for a given value of the conduction 

parameter, the resulting effectiveness will not be affected by small 

* 

changes in As . 

From the schematic representation shown in Fig. lb, it should be 
noted that the left edge of the regenerator is physically the same os 
the right edge and, therefore, the matrix inlet temperature for the 
elements of the first colutan on the side of is the same as the 

outlet temperature of the corresponding element of the last column on 
the side of Expressed mathematically. 



where i. and f can take the values of l, 2, 3, b' r . 

This is referred to as the reversal condition. 

From equations (3), (4), (5), (9), (10), and (11) It is seen that 
in order to solve for the outlet matrix temperature for a particular 
element, it is necessary to know the matrix temperatures of the next 
element. An estimate of these values can be made by determining the 
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temperature distribution for the ease of zero conduction in the direction 
of fluid flow, since the constants D 3 and Dg in equations (3), (4) and 
(5) and the constants F 3 and Fg in equations (9), (10) and (11) are zero 
in this case. 

A temperature scale can be used for which the fluid entrance tempera 
ture is zero at Che side of C n j n and unity at the side of C taax . 



r ‘(i,g) " ° 


for g = 1,2,3, 


K n 


x(l,j) - 1 


for j = 1,2,3, 





For calculation purposes a temperature distribution is assumed on 
the left edge and the problem is solved for the no conduction case in 
order to obtain an initial estimate of matrix temperatures for working 
the problem with conduction. This also provides the no conduction 
effectiveness necessary for comparison. 

The outlet temperatures are calculated for every element by repetitiv 
use of equations (3), (A), (5) and (&) depending on the Location of the 
element. The calculation is started with the first element, first 
column on the side of C lLax and then working down the column. When the 
first column of elements is completed, the second column is calculated 
and so on until the outlet temperatures of all the elements on the side 
of Cmax are determined. It was found that convergence could be enhanced 
if two passes per column were made before proceeding to the next. one. 

The reason for this is that the equations for the outlet temperatures 
contain the matrix temperatures of the next element which are necessary 
to estimate. 

At the seal represented by the double line in Fig. lb, 



T x(i,N x + 1) - T n ( f , l ) 



since the matrix outlet temperature for the elements of the last column 
on the side of C ma y is physically the same as the inlet temperature of 
the corresponding element of the first column on the side of C^ n . 

The 3arae method is applied to the side of as to the side of 

Ctsjax* If the temperature distribution assumed on the left edge was 
correct it would then be duplicated on the right (i.e., the reversal 
condition is fulfilled). If, however, this is not the case, the 
resulting temperatures on the right side .T n (f f u n are now used on 

the left T x (i,l) and the procedure is repeated, After each pass (i.e., 
the complete calculation of a temperature distribution) an energy 
balance is made; and before the solution is accepted for a particular 
set of parameters, the heat balance error together with the reversal 
condition has to be fulfilled to the specified accuracy. 

The effectiveness of a heat exchanger is defined (Appendix 2) as 



From the conditions of the problem the above expression for effectiveness 
reduces to 





The heat balance error was computed from 



1 . v AV 1- JO £ 
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That t’ 5t>).4it i ’.'ill ••O'lVC.Ti'C (I.-,, 



«. i<a revere l conditio:, n.t 



and this heat balance *rror converge to zero) will depend on the number 
o £ elements used. More Important than this is the relation existing 
between the numb.v of subdivisions for the three streams. The sufficient 
but not necessary conditions for convergence are (Appendix 1) that 



N\ 



\ A * ' 



\\ 



Y 
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Fig. 2 is a plot of the different values of effectivencoo obtained 
for the following set of parameters- 

C . /C - 1.0, Cr/CL,„ » 2.0, (hA)* - 1.0, NTUo - 10. 0 \=0.1, As" = 1. 

Several curves are obtained depending on the relation existing between 
the number of subdivisions for every stream. To e:ctrapolate the 
effectiveness to an infinite number of elements, these curves nay be 

f 

approximated by straight lines. The error introduced will depend on 
the number of elements used and will not be constant for all the range 
of the parameters. As on illustration consider the set of parameters 
used to plot Fig. 2. Three cases arc investigated: 

1) tJ :? = ir n - N r /2 

2) •?- - \n - Uc 

3) V. - n u - 2H r 

A linear extrapolation is performed using the values of effectiveness 



obtained for values of 



Froo the results it is found 



N r (N x + N n ) 

that the minimum N r necessary to give an error in the cxtropolatcd 
effectiveness of less than 0.01Z is 18, for which » 



76.57% 




(\J 

to 

•H 



1 A 



l/N r (N x+ N n ) 




for the three cns< s. The determining faetor in the extrapolation is 
the nunber of matrix elements. The higher the value of '. r U 3 ; d for the 
extrapolation the smaller the error in the straight line approximation. 
Unfortunately, an increase in M r will bring a higher increase of K x and 
N n , since the convergence condition for this example is 

N 

\» * 

That this is only a sufficient condition can be seen from the case of 

Nx - = Nr/2 

for which the solution converges for values of N r >5; and when 

»x " N n - N r 

the solution will still converge for N r >20. It is true, however, that 
the closer the number of matrix elements are to the sufficient convergence 
condition, the fewer passes it takes for the solution to converge to 
the specified error in the effectiveness. Table 1 shove the number 
of passes necessary to obtain a . olution with a heat balance error of 
less than 0.0057,. 

For a given set of parameters the effectiveness vas determined 
first using one set of subdivisions and then recalculated increasing 
the number of subdivisions. A linear extrapolation was performed to obtain 
the effectiveness corresponding to an infinite number of elements. The 
accuracy for each value of effectiveness was specified to be of four 
significant figures. When these values were extrapolated to an infinite 
number of elements, the accuracy vas reduced to only three significant 
figures since a compromise had to be made with the time available for 
computations and th<.: number of runs to be made. However, the method 






TAolA l 

Number of passes necessary to obtain an effectiveness with a heat 
balance error of less than 0.005* for the following at t of parameters : 
<«ln/Cpax -1.0 . Cj/Cn^, - 2.0 , <hA)* = 1.0 , 



HTUo - 10.0 , X* °‘ l » As * -1.0 



Hr 


V 


*»n 


Mo. of passes 


vV 

7, 


12 


12 


12 


27 


76.581 


12 


26 


24 


22 


76.581 


14 


7 


7 


41 


76.578 


14 


14 


14 


36 


76.576 


14 


28 


28 


27 


76.576 


16 


16 


16 


41 


76.574 


16 


32 


32 


32 


76.574 


18 


9 


9 


62 


76.573 


1.8 


18 


13 


50 


76.572 


10 


36 


36 


38 


76.572 


20 


20 


20 


60 


76.572 


j 22 


22 


22 


71 


76.571 



* Ojo is obtained by a linear extrapolation of the effectiveness 
for N c and N r -1. 
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imposes no restriction as to the ultimate accuracy of the solution. 

When working with very high values of C r /C, n ^ n the computation time 
becomes considerably long. This is due to the fact that the final 
temperatures differ very little from the assumed initial temperatures and 
a large number of elements is required. The starting estimates of matrix 
temperatures then become very important. For the limiting case of C r /Cjj,i n 
equal to infinite, there is no temperature variation in the direction 
of metal flow so that the method without modifications is not applicable. 
For this case, however , the exact solution of Hahnemann 3 is available. 

The Fortran system was used to program the above method for the 
CDC 1604 digital computer. The matrix notation employed in the equations 
derived previously is very convenient in writing up the program, since 
several loops can be made by just changing the subscripts of the variables 
involved. The complete computer program is given in Appendix 2. 



1 ? 



3. Discussion o£ results 



The solutions generated ace presented in tabular form in Tables 2 
through 25 ancl graphically in Figs. 3 and 4. The accuracy in the 
effectiveness was specified to be of three significant figures. It was 
found that for 

0.7 (hA)*^ As* ^1.3 (hA)* 

the resulting values of effectiveness were almost identical to the ones 
calculated for 

As* = (hA)* 

so that these were the values used in preparing the Tables. 

By using N’fUo and X os defined in the nomenclature, which contain 

(hA) and As' explicitly, the influence of these two parameters on the 

• -A C? 

conduction effect, , were found to be very small for values of 
Cjj^n/Ctuax "> 0.9 and \^0.1, which is the range for the gas turbine 
regenerator problem. 

Figs. 3 and 4 show the effect of conduction on effectiveness for 
values of Cmin/Cmax = 1.0 and 0.9 respectively, with different values of 
NTUo, and \ , up to a limit of 10%. The conduction effect 

continues to increase with increasing X • For small values of A this 
effect increases nearly exponentially. By further increasing X » a 
saturation is reached which approaches the limiting condition of \ = oo. 

Fig. 5 is a representative curve whLch shows the limiting influence on 
the conduction effect by continuously increasing X for a given set of 
parameters . 

The conduction effect for values of C r /C^ n = 10. 0 arc almost identical 
with the values reported by Schultz 2 and Hahnemann [3 ] , 6 for the 
limiting case of = oo. 

io 



l'or email values of NTUo the effect of increasing Cr/Cciin ^ to 
increase slightly the conduction effect, as can be seen from Fig. 6. 

For higher values of NTUo the effect of conduction is reduced with 
increasing Cr/C,^. 

The effect of increasing NTUo upon is intimately related to the 
value of Cp/Cnin and C in £ n /C raax . This effect can best be explained by 
observing Fig. 6. For the case where C^n/C^ax =1.0 the effect of 
increasing NTUo will result in an increase in 4r regardless of the value 
of C r /C, n ^ n ; of course, the higher the value of C r /C m ^ n the smaller this 
effect will be for relatively high values of NTUo. Now, by reducing 
^min/Snax sn< * increasing C r /G [|t ^ n a point is reached in which an increase 
in the NTUo will result in a reduction in ^ . The point where this 
occurs will depend upon the value of the conduction parameter X . 
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A. 



Surnr *& r y * : ad Co ncl us ions 

A finite difference numerical analysis is presented which v/ill 
determine the effect of the longitudinal heat conduction in rotary 
regenerators for steady state conditions. 

The ranges of governing parameters which have been covered are; 

0.5 ^Cfoin/^max 
l.O^C^^^JO.O 
1.0 6 iJTUo ^20.0 
0.25 4(hA)*4l.O 
0 4 A =, 0.2 
0.15 4As*41.3 

The ranges of parameters of interest for gas turbine regenerator 
application are plotted in Figs. 3 and 4. A more accurate determination 
of the conduction effect than that given in the figures can be made by 
interpolating values from the Tables. 

The computer program used to carry out the calculations on the 
GDC 1604 digital computer is included in Appendix 2. This program 
imposes no restrictions on the parameters for which the analysts would be 
applicable. It is recommended, however, that for values of Cr/C m ^ n >10 
the solution for C r /C m £ n - 10 be used. For this case the computer time 
for a solution increases considerably and the error introduced by assuming 
Or /Onin = 10 is negligible. 
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EFFECT 
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- 1.0 
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. * 
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1.0 






e % 


■ 




Conduction effect,— 


7. 






NTOo A 


0 


0.01 


0.02 


0.04 


0.08 


0.14 


0.16 


0.20 


0.32 


l 


46. 64 


0.42 


0. 78 


1.48 


2.61 


3.94 


4.30 


4.95 


6.42 
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60.06 


0.50 


1.11 
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3.81 


5.93 


6.54 


7.65 
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3 


66.72 


0.68 


1.32 


2.50 


4.58 


7.17 


7.92 


9.30 
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4 
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0.77 


1.48 


2.82 


5.17 
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8.93 
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5 
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1.63 


3.09 


5.65 
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15.44 


6 
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6.05 


9.42 
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7 
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1.87 




6.41 


9.93 
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8 
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| 


1.97 




6.72 


10.33 




13.34 




9 
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2.07 




7.01 
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7.26 
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14.24 
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12 
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7.71 
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14 
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16 
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18 
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20 
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40 
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3.55 
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i 
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4.78 


7.97 
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1 
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4.04 
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12.67 


16.13 


13 


89.37 

i 


2.54 


8.44 


12.83 


16.31 
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2.59 
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6.96 
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7.59 
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15.12 


90.18 
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2.77 


4.16 
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i 


18 
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1 
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1 
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20 
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1.81 
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TABUS 9 LONG ITUD1 HAL 


HEAT CONDUCTION EFFECT 






^In^Snox “ l -° 


» Cf/Prain 


- 10.0 , 


(hA)* * Aa* 


- .5 






e % 


Conduction ef£ect ; ~ % 


NTUo 


0 


.015 


.03 


.06 


.12 


.24 


.48 


1 


49.96 


0.65 


1.23 


2.24 


3.84 


6.04 


8.51 


2 


66.59 


0.90 


1.72 


3.20 


5.67 


9.32 


13.84 


3 


74.91 


1.02 


1.97 


3.70 


6.62 


11.04 


16.64 


4 
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APPENDIX I. 



Derivation of iterative equations. - 

Equations (1) ami (2) from page c can be put in the following form: 

t — f -Cl . T — T — T (l.l 

*(•« x ' s -» + j 1 '-'x’.jj 'H '-* 1 .}) *v.<. j; V. j*') 



" v<-ti y, 




V- ^ ) 







( 1 . 2 ) 



where 




N I 




J 



I 



/ \ 




and 

V— * 

Equations (1.1) and (1.2) may be written as 







(1.3) 



and 



r U 



(1.4) 
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i)y a tiding to equation (1.3) the product of ( i+C^ ) ( l . 4) and 



(*) 



where 



and 

- J-i 



From equation (1.1) 



( 6 ) 



where 



ft 



and 



For the elements from the first and last row equation (1.2) becomes 



and 



(1.5) 



( 1 . 6 ) 

\ . - i 
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respectively. T w t *e same procedure ns before mid solving for the matrix 
outlet temperatures, the following equations are obtained: 




-t 



■). i 



-4 



(3) 



and 



h 



(5) 



for the first and Inst row eLeraents respectively, where 



and \J - 

f 



Equations (7) and (8) ntay be written as 



(1.7) 



and 




( 1 . 8 ) 



Equation (1.8) reduces to 
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t -i 
^w-, S ) Wi, a ) 

and 



U ‘3 / 'd' 



+ t 3 



’ '-n _‘y^ _-p 

(^’ ;}) ^ U ^r,3*H) ° i ^t*“ , 5) ^W r “ ,g+l)j (1.9) 






- I , 



( 1 . 10 ) 



for the first and last row elements respectively, where 



- N I 



- V 



By comparison of equations (1.7) and (1.8) with equations (1.1) and 
(1.2) it follows that the form is the same and only the coefficients are 
different. The same holds true for equations (1.9) and (1.10) with 
equations (1.6) and (1.5) respectively. By substitution the following 
equations are obtained: 



o . .. V; 



’ < h itf. i) 



(9) 



'.Nr’’ 



i : 



I. 'I. 



( 10 ) 



I • 1 



, 1 






(U) 



for the first, middle and last row elements respectively and 
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( 12 ) 



which is the same for oil the elements, where 



I 



i 



; 



and 
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the relieving form: 



Conv^egsucc eouJitior, - 

Lquation (i *2) car b~ in 



f 



(i.ii) 

n tha energy balance equations, the aide ot Crn 3x was assumed to be the 
hot side so that by the second law of thermodynamics 



and from equation (1.11) 

t- I . — 



Rear?.*, tv i *r equation (1.12) 



( 112 ) 



(1.13) 

Again by the second law oi thermodynamics the matrix temperature variation 
in the direction of fluid flow is ouch that 



(l.U) 

<( ^ 
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aUd 



(1.15) 



so that 



(1.16) 



Substituting equation (1.16) in equation (1.13) 



(1.17) 



Three cases are investigated: 

1. - If C 3 C, then equation (1.17) becomes 

> — 

j * 

2. - If Cj » 0 then equation (1.17) becomes 

y j ■ • 

3. - If C, • 2C^ * 0 then equation (1.17) becomes 



From the above cases it can be seen that a sufficient condition for 
equation (2) to converge is that 



Substituting the values for the constants and rearranging yields 
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For the first and last row elements the convergence criterion 



becomes 



which for substitution of the parameters gives 



This condition is less stringent than for the middle row elements. 

Similarly, Cor the side o. Crain the sufficient condition 
for the middle row elements becomes 



which for substitution of the parameters gives 



For the first and last row elements it is found that 



which requires that 



is the sufficient condition. This again is a less stringent condition 
than that imposed by the middle row elements. 
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r» i i (even t ial equ.it io u . - 



Hfius* n 4 hi*s given a theory for regenerator* without heat 
conduction. When the wiit ussumpt ivr.tt -re nr. Se, but tin. longitudinal 
conduction of heat is included, (See page * ) the differential 
equations obtained are _ : 

for the aide of C max (gas flowing from ; - 0 to * L) ; 

i 



and for the side of 'reversed gas flow) 



A temperature scale can be used for which the gas entrance 
temperature at the side ol C T1 j n is zero and that at the side of C 

MS 

is unity. The boundary coi.a tions are then 



and 



APPENDIX 2 



Computer program. - 

The Fortran system was used to program the problem for the CDC IbO'f 
digital computer. 

The sij: dimensionless parameters are used as input together with 
increments and factors to change the values of NTUo, and As, so that 
solutions can be generated for several combinations of the parameters 
without having to feed into the computer these different sets of para- 
meters every time. 

Another input is the number of subdivisions of the three streams 
together with the increments, so that when the computer finishes the 
first series of runs for a given number of subdivisions it automatically 
increases the subdivisions by £he specified increments and the two values 
of effectiveness obtained are extrapolated to an infinite number of 
subdivisions. 

In order to obtain the initial estimates of matrix temperatures, 
the problem is first solved neglecting the longitudinal heat conduction 
as explained previously. For this case the equations derived Ln Appendix 1 
are simplified to (for detailed derivation see Reference 1) 
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The inlet matrix temperatures necessary to initiate the 

problem for zero heat conduction are estimated by assuming \, , 

(See 1'ig.l'. ) rnd feeding this value into the computer. The retraining 
inlet matrix temperatures are obtaLned from the relation 

H ‘ ij ~ 'i ' 

The number of pusses necessary to meet the desired accuracy in the 
solution changes accordingly with the initial estimates of matrix inlet 
temperatures, but considering that the computer time required to work the 
problem for zero heat conduction is so small as compared to the case with 
conduction, that these estimates do not become important. Only for high 
values of Ci/Qain these estimates may be of significance since for this 
case the change in the reversal condition after every pass is very small. 

In general, it was found that an estimate' T, jv 



of 


0.4 to 0.6 


for 


^r /^in 


<5.0 


and 


0.7 to 0.9 


for 


Cr/Cmin 


^ 5.0 



together with the above relation for T x ,^ should converge to the 
solution in a reasonable time. 

The temperature distribution is calculated by repetitive use of the 
iteration equations. 

The effectiveness may be defined as 



0 Hi ax ( ^ XoJ ~trio ~ 

C-n viw (t^xi — T-n,*) Txi - Vi 
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irjio ' s- t j oiv; of i> 



-rcMcu 




therefore, the expression for effectiveness for the definition reduces to 



N , 




which is the expression used for the computation of the effectiveness. 
The heat balance error is computed from the relation 

( A , - u ) 



error = 



J. V 






which reduces to 



error = 
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The extrapolated of Cectiveneas is conput«d fren the relation 



(iv 



2 ) 



i 



^ i 



where 



The input parameters tn Fortran notation are 
^ = Sain^aax > i 2 - C r /C, n j n , F3 = (liA) , 

PA = NTUo , P5 « \ , ?6 = As* 

The initial values of NTUo, \ and As* are changed by the following 
relations 

?A « pA PAIN 
P5 = P5'*?5FA + P51N 
P6 = P6 - p6lN 

where PAIN, P5FA, P5IH, P6IN are also input values. The inputs pAFI, 

P5FT. and P6FI are the final values of PA, P5, and P6 respectively. 

The program is written so that the computer after working the 
initial set of parameters, i' reduces '6 by P6IN every time until 
P6 becomes less than P6FI. Then it sets F6 equal to the initial value 
and changes P5. The procedure is again repeated until P5 becomes greater 
than P5F1. Then P5 is set equal to the initial value and PA is increased 
by PAIN. The cycle is repeated until PA becomes equal to or greater than 
PAFI. 
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I 1 INPUT VALUES J 
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RE • VARIABLE USED TO EXTRAPOLATE EfTECTI VENESi 


EM2=NM2 
1)06 L= 1.NM2 
t -t 

RT*1.0-t(E— 1.01/EM2) 
6 UH2C L 1 =UH ( I, Il'KT 


ESTIMATES Or T x ( lf l) TOR 2nd RUN 


7 FORMAT!// I4HN0 OF E L E MEN l S3X , 3 l 5 , 4 X , 3 ( 5 ) 

WRITE OUTPUT TAPE 2, 7 »NM,NH»NC,NM2, NH2 * NC2 
PRINT 7,NM,NH,NC,NM2*NH2,NC2 




8 I C= 1 


INDICATES lit RUN, USING INITIAL No Or 3UBDIV. 


9 EM=NM 

EH = NFf 
6C«NC 


riXED TO FLOATINO POINT CONVERSION 


NMM=NM- ( 
NHP’NHF ( 
NCP=NC+I 

NM P = NM ♦ I 


t 

DEFINITION OF NEW VARIABLES 


00 10 J* 1 ,NH 
10 TH ( 1 * J 1 = 1 .0 

UO 11 L-1.NC 

1) T C C 1 , L 1 -0. 0 . . _ 


SET > 1.0 AND t n (l,g) * 0.0 


) 2 A 1 0 = (.0*EM/(PI«P2»EH1*(2.0 - EM*P31/(P4#PI«( 1 . 0* P 3 1 ) 

A 1 = 2. 0/A 10 

A20= 1.0*1 p 1 • P2 • E H 1 /EM+ (2 ,0«EH»P3*P2 l/(P4.(|.0*P31l 
A2=2.C/A20 

A30=1.Q*EM/(P2»ECI ♦ ( 2 .0 • E M ) / l P4 • < 1 . 0 ♦ P3 1 1 
A3-2.0/A30 

A40= 1 .0*(P2»EC 1/EM +(2.0*EC«P21/(P4»( 1 • 0* P3 1 ) 

A4=2 . 0/ A40 


CONSTANTS FOR THE NO CONDUCTION CASE 


C HOT SIOE NO CONDUCTION 




121 M a 1 


M DENOTES PASS NUMBER 


13 00 ) 5 J=I.NH 

DO IS 1 = 1* NM 

THClF(,Jl = THCI t JUAI-CUH(l,Jl-IH(|,Jll 
15 UH ( I * J * 1 1 = UH ( I , J 1-A2» CUHCI ,J|-TH(I , Jl ) 


ITERATION ON SIDE Or C B *x 


C COLO SIDE NO CONOUCTION 

00 17 K= 1 * NM 

17 UC(K , 1 1=UH(NM+ 1-K*NH* 11 

UO 19 L = 1 » NC 
0019 K=(,NM 

TCCK* 1,L) = TC(K,L 1*A3» CUC (K.U-TC (K.L 1 ) 
19 UC ( K , L ♦ 1 1 =UC ( K , L I-A4MUC 1K,L)-TC(K,1>) 


ITERATION ON SIDE OF C»i D 


01=0.0 

00 20 J= 1 * NH 
20 GI>TH(NM«1.J1«GI 
THO=Gl/tH 


TH0 • [tz(K r d,J)j Brg 



G2=0 .0 
00 2 1 L = 1 1 NC 
G2*TC( NM* 1 ,L I + G2 
TC0IMHG2/EC 



TCC - EFFECTIVENESS 



ERR(M1=1.0-IP1«TCO!MI )/( 1.0 Til I 
IFIIRR(M) )23,24,?4 



COMPUTE HEAT BALANCE ERROR AND COMPARE IT 



c««ini*-cnn»ni 

24 IE(ERR(M)-4*0f-6l 30, 30 ,25 


WITH MAXIMUM SPECIFIED 


1 2S (FI90-MI 26,26,215 




j 26 IF 1 ERR! M- 1 1-ERRI M 1) 2 7,27,2? 




| 27 PAUSE 1 1 


CHECK CONVERGENCE 


1 28 IFI SENSE SWITCH 1 1 32, 121 


PAUSE 11 = SOLUTION CONVERGES 


27 PAUSE 12 


PAUSE 12 = SOLUTION DIVERGES 


GO TO 28 




215 DO 22 1=1, NM 




22 UH 1 I , 1 1=UC 1NM+ 1-ItNC+l ) 


REVERSAL CONDITION 


M-Ht 1 




GO TO 13 





301 

302 



T 2 = TCO l M- I )- ICO ( M J 
1FIT2) 301,302,302 
T2= - I 2 

IF ( T2-4.0E-61 305, 305, 2 ‘j 



305 

31 

312 



1 F I (C- I 1 31*31, 33 
DO 312 J=1.NM 
UH 1 I J )=UH( J , 1 1 
6 F C 1 , 1 1 = TCOl M 1 
ER * 1, I )=ERRt M 1 
GO TO 3V 



CHECK DIFFERENCE IN EFFECTIVENESS IN 
THE LAST TWC PASSES 



STORE T x ( itl ) , EFFECTIVENESS AND 
ERROR FOR THE 1st RUN 



FORMAT ( 4FB.2, 1 5, F 10.6, F7.2, » 31 

PRINT 32,PI,P2,P3,P4,M,TCO(*>,ERklM) ,!C 

P^USE 13 

IFCSENSE SWITCH 1) 32 1 , 1 2 1 

IF 1 SENSE SWITCH 2) 1,322 

NM=NM2 

NH=NH2 

NC=NC2 

G0_T0 _5.. 

EFF=TC0|M1* 

B1=f FF 



INCREASE NO. or SUBDIVISIONS AND REPEAT 
1st RUN 



l.O+RE 1-EFI I , I l«Kh 



1 34 


FORMAT! //4F8. 2, 




WRIIE OUTPUT T l 




00 35 L= 1 , NM 


35 


UH2 I L 1 =UH1 L , 1 ) 


39 


I A=2 




I B= 1 



+ 9. 21 ,F 10 .6) 

/ ,P 3 , P 4 ,EF I 1 i 



EXTRAPOLATE EFFECTIVENESS AND STORE 



1) , E K ( 1,1) , Tl 0!*) , E*R(*> , 



STORE T x ( it i) FROM 2nd RUN 



40 DO 44 1=1, NM 

Ml 00 42 J- 1.NHP 

42 UHC { I , J 1=UH( I , J ) 

43 00 44 L=I,NCP 

44 UCC! 1 ,1 HUCI 1,L 1 

C CONSTANTS 

50 £1=(P4*1 1.0+P3) 1/12. 0*EM ) 

61 1* I.0+E1 

F 2= C P2«EC)/EM 

E 3=0 . S*PS*EM 

E4=E 1+{ E 1 1*1 2.0*E3+E21 ) 

E5 = E HIE 1 I •( E2 + E 3) 1 
F 1=1 2.0*6 1I/F4 

F2=( E 1 + 1 El 1*1 2.0*E3-E2 )) 1/EM 

F 3* 1 E 3* E 1 11/EM 

F4 = I 1.0-E I 1 / 1 1 1 

F5 = E I/El 1 

F6=( 2.0*E 1 1 /t 5 

F7=JE 1 ♦( E 1 I • ( E 3-E 2 ) 1 |/kS 

F8= ( E 3*E 1H/P5 

C1=1P1«P4*1 1.0+P31 )/l?.0*P3*tM) 

Cll=1.0+C1 

C2=1P1*P2«FH)/E M 

51 C3-1PI*PS*EMI/I2.0*P61 
C4*C1 + (CU«(C2 + 2.0*f 3) > 

CS=C !♦( C 1 1 * I C2+C 3 1 I 



SUBSCRIPTS FOR ERROR AND EFFECTIVENESS 



UHC(I.J) = T x ( i , J } , UCC(I,L) . Tq{ f f g ) 
FOR CONDUCTION CASE 



CONSTANTS FOR CONDUCTICN CASE 



5 3 
C 

60 

6 ) 

62 



64 

66 



67 

60 

69 
C 

70 

71 

72 



73 

74 



75 



D1*(2. 0*C1 )/C4 

D2=IC1*(C1 Nt2.0*C3-C2) ) J/C4 

03* I C 3*C 1 I )/C4 

04= ( I . 0-C ll/CII 

05 = C 1 /C I 1 

06= I 2 «0*C I 1/C5 

07= I C I ♦ I C 1 1 *1 C3-C2 1 )J/C5 

08= I C 3*C 1 1 1/C5 | 

M= I 

HOT SIDE WITH CONDUC 1 I 0»M | 

00 67 J»1,NH 
L= ) 

UHCt l , J+ 1 ) =1)6 • T H I 1 , J )-0 7*UHC I 1 , J ) *00 * ( UHC I ? , J 1 ♦ UHC ( 2 , J * I 1 » 

TH( 2, J )=04*THt I, J ) *05* ( UHC ( I , J 1 *UHCI I , J* Ml 
0064 1=2, NMM 

UHC ( 1 , J*1)=U)*THt 1 , JI-02HJHCI 1 , J I *03 • ( UHC ( 1- 1 , J ) *UHC( )-) , J* 1 1* . 

I UHC ( 1 * 1 , J 1 *UHC I IH.Jtlll 
TH( ) * I , J 1=04*TH( I , J ) *05* l UHC t I , JJ+UHCI I , J* III 

UHCINM, J* 1 )=06*TH( NM, J )-C7*LHCINM,J) ♦ D8 • I UHC l NMM , J ) *UHC( N*»m, J* 1 ) ) i 
THl NMP, J»=04*THlNM, J J +D5MUHCINM, J )*UHCINM, J* 1 J ) 

L = l* 1 

) F ( 3-L ) 67,67,61 

CONTINUE 

DO 69 K=),NM 

UCC1K, I 1 =UHC ( NMP-K , NHP ) 

COLO SIDE WITH CONDUCTION 

00 75 L=1,NC 

J=1 

UCCI 1,1*1 1 =F 6*TC I 1 , L ) - F 7 »UCC 1 1 ,U *FH • ( UCC I 2 , L ) *UCC(2 , L + 1 ) 1 
TCI 2,l)=F4*TC< 1 , L J ♦ F5 * ( UCC I 1,L)*UCCI 1,1*111 

00 73 K = 2 » NMM 

UCCI K , L ♦ 1 I =F 1*TC(K,L)-F2* UCC I K , L 1 ♦ F 3 * I UCC I K- I ,L> *UCC(K-I ,L* 1 )* 

1UCC IK* I , L 1 *UCC I K* I ,L + I 1 1 
TCIK* 1 , L ) = F4*TCI K, L 1 * F 5* I UCC I K , L 1 ♦ UC C I K , L * I ) 1 

UCCI NM,L* l J = F6*TCI NM,L J-F 7 »UCC I NM , L I *f 8* I UCC I NMM, L 1 *UCCC NMM, i * 1 1 1 
TCtNMP.L }=F4»TCI NM , L ) + F5 • I UCC I NM , L 1 + UCC I NM , l ♦ I) 1 
J=J* I 

1 F ( 3- J 1 75,75,71 
CONTINUE 

G 1 =0 « 0 



DO 76 J= 1 , NH 
76 G1=THI NMP, J 1*G 1 

TH0 = G1/EH 

G2 = 0 • 0 

00 70 L = 1 , NC 



70 G2=TCINMP,L ) *G2 

TCP I M 1 =G2/FC 

ERR { M I = 1 .0*- 1 P 1 * T CO l M J )/( I.O-THOI 
1 F < ERR I M ) ) 00,00,01 
00 ERR 1 M ) = -ERR | M 1 
81 1 F I ERR I M )-4 . OE-6 ) 90, 90, 6? 

02 IFI90-M) 03,03,05 
83 IFtERRIM- 1 >-ERR<M) 1 04,84,87 
04 PAUSE 22 

04 1 I F I SENSE SWITCH II 042,53 
042 1FISENSE SWITCH 2) 00,084 
85 M=M* 1 



00 06 1=1, NM 

06 UHC I I , I 1 =UCC l NMP- I , NCP I 

GO TO 60 

07 PAUSfc 21 
GO TO 04 I 

80 1 F I IC- 11 001,001,082 

00 I EF| I A , 181=0,0 
ER( I A, 101=0.0 
GO TO 97 

002 FORMAT | 1 1HS0L OIV FOR I 3 , I 0x , E 9 . 2 , F 0. 2 , 2 I F 1 0 . 6 , F 9 . ? 1 1 

WRITE OUTPUT 1 APE 2 , 882 , IC , 05, P6 , F Eli A , I 8 1 , FR I 1 A , l f! 1 , TCO I M ) , 
I QO TO. 97 



RR ( M 1 



CONSTANTS FOR CONDUCTION CASE 

_ J 



iteration on SIDE OF C..X 
TWO PASSES PER COLUMN 



ITERATION ON SIDE Or Cain 
TWO PASSES PER COLUMN 



884 FORMAT I 4F 8 ♦ 2, E 8. 2, F 8.2,2*, |$ V FI0.6,E9.?1 

PRINT 884, P I,P2.p3,P4,PS,P6, 1C.TC0IM ] , F RR ( M 1 
PAUSE 23 

IFISENSE SWITCH 1) BPS. S3 
BBS IFISENSE SWITCH 2) 114, 8e 

89 IFIIC-11 891,891,93 

8 9 I EFUA. I 8 I = TC 0 1 « 1 

ER I |A, I B 1 =ERR ( M 1 

GO TO 102 

90 TI*TCOIM-l 1-TC0IMI 
I F | T II 9 1 ,92,92 



91 T 1 — — T 1 

92 IF ! T 1-4 • OE-6 I 89,89,82 



EFF = TCO(MI*( I.O*REI-FFIH,IM»Rb 
82* C I B 1-6FF ) • 100.01/P 1 



97 

971 

972 

973 

974 



FORMAT 1 32 X, E8.2.FB.2.2IF I0.6,E9.?I,F10.6,F8.31 

WRITE OUTPUT TAPE 2, 94 , PS , P6 , E f 11 A , 1 8 1 , fc R I I A , I 8 1 , TCO I M ) , E RR l M I , fc FF 
1,82 

IF I SENSE SWITCH 3) 9S, 102 

FORMAT I F 8 . 2 , E9.2.FB.4 1 

PRINT 95 , P4, P5 , 82 

GO TO 102 

00 974 I=),NM 

00 972 J= I , NHP 

UHC ( 1 , J )=UH{ I , J ) 

DO 974 L- I , NCP 
UCCt 1 ,U=UC( I , L ) 



102 IFISENSE SWITCH 3) 131,132 

131 FORMAT! /4F8.2,E8.2,F8.2, 3x , I 3/ ! 1 SF 7. 4 1 ) 

WRITE OUTPUT TAPE 3 , 1 3 1 , P 1 , P 2 , P 3 , P4, PS , , I C , 

UUHCt J, 1 1 , UCC ( NMP— J ,NCP) ,J = 1 ,NM 1 

132 P6=P6-P61N 



18= IB* 1 

1 F { P6-P6F 1 ) I03.103.S1 

103 P6- S6 
I B= I 

104 PS*PS*PSFA*P5IN 



1 A= 1 A ♦ ) 

1FCP5FI-PSI 10S,S0,S0 

105 P5-S5 

1FIIC-1) 106,106,108 

106 IC=2 



NMFNM2 

NH*NH2 

NC*NC2 

DO 107 K=1,NM 

107 UHIK, I )=UH2IK) 

GO TO 9 

108 I F 1 P4F I— P4 1 113,113,109 

109 P4=P4*P4 IN 



NM-NM-NM 1 
NH* NH-NH 1 

NC = NC-NC 1 
00 111 K- I » NM 
111 UH1 K, I I-UHIIK ) 

__ GO TO 6 

113 PAUSE SO 

IFISENSE SWITCH 11 114,1 

114 DO 1 IS L * I , NM 

115 UH ( L , I 1 -UH2 1 L 1 
GO TO 5 

ENO 

ENO 



ro 




UNCLASSIFIED 
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